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ABSTRACT 



Loran-C time difference (TD) readings were taken 
consecutively at several survey stations near Monterey, CA, 
over a period of three days. One station, Range-7, was 
designated the "known" point of a differential Loran-C 
system. Readings from the known point were used to correct 
readings from a second survey station. Luces Point. A method 
of improving the precision of Loran-C TD readings based on 
the redundancy and relative accuracy of three LOP's was 
developed and applied to the data. Since only one receiver 
was available, a linear regression of TD vs time was 
calculated and used in the differential correction. 

Based on 496 sets of data taken at 5-second intervals at 
Range-7 in two groups, before and after about 250 readings at 
Luces Point, the absolute accuracy of Luces Point data was 
improved from about 385 m to about 48 m compared to the known 
position of the point. Precision was improved from about 
14.9 m to about 12.6 m circular error (CEP) using the three- 
station correction. Further improvement would probably have 
resulted if two receivers were available for real time 
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i. DfrooDUcncN 



A. Purpose 

The purpose of this study was to investigate the potential 
accuracy of the Loran-C navigation system using a differential 
correction technique and a correction for the use of three secondary 
stations rather than the usual two. 

Sources of error were first examined and their magnitudes 
estimated. Methods of expressing the accuracy of a determination of 
geographic position were also considered. 

It was determined that position accuracy is composed of two 
factors: the absolute accuracy of the position (which can be 
considered the bias, or systematic error, in the system) and the 
precision of the individual readings. 

Over longer periods of data collection it was considered possible 
to subtract out the system bias by means of a quasi-differential 
correction based on intermittent time difference readings at a fixed 
point, designated the "known" point of the differential system, and 
to improve the precision of TD readings by taking advantage of the 
redundancy of the reception of three master-secondary pairs. 

For this purpose, a theory of corrections was developed and an 
experiment was designed to determine the effectiveness of the 
corrections. Hie experimental data was processed by computer, and the 
relative accuracies were calculated. Computer programs for processing 
and displaying data were written for this purpose. 
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B. Method 



To collect data for processing, several known survey points with 
published positions were chosen as locations for collecting Loran-C 
data. Using known points provided a means of comparing results. One 
point in a relatively central location was chosen as the "known" point, 
where a larger amount of data was collected. Data from the other points 
was corrected based on the data from the "known" point. 

During periods of about eight hours, the receiver was moved from 
one survey point to another. The antenna was set up over the points 
and time difference data was collected for periods of about 20 to 30 
minutes at each point. The time required to move the receiver from one 
point to another varied from about 30 to 60 minutes. 

Data was collected by a Racal Megapulse monitor receiver with 
internal clock accuracy of about 0.02 fj.s ec. A Silent 700 data terminal 
was used to collect the data on cassette tapes. Data was recorded at 
intervals of 5 seconds. Later, the data was transferred to 5 V disks 
for use on a computer. 

Since there was no way to collect time difference data at two 
points simultaneously, a method of estimating the time differences at 
the known point was developed. This involved using the available data 
to establish linear regression lines of time difference versus time at 
the known point. 

Using the known time of data collection at the other points, an 
estimate of the time differences at the known point could be used to 
calculate a correction for data at the other points at the times of 
data collection. 
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A number of computer programs were written to calculate the 
required information for differential corrections. A Correction For 
Three Secondary Stations was developed mathematically and used to 
process the data in an effort to reduce the random error of the time 
difference readings. 

Based on reduction of standard deviations of the data and 
reduced circular error, the three station correction was 
effective to some degree, although absolute error turned 
out to be the limiting factor in the experiment. Reduction 
of absolute error by means of the differential correction 
was significant. 
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II. THE LORAN-C NAVIGATION SYSTEM 



A. Basic Description 

Loran-C is an electronic positioning system based on the reception 
of pulsed LF (90-110 kHz) signals from a master transmitting station 
and one or more secondary stations which comprise a Loran-C "chain". A 
brief description of the system will be provided here as a basis for 
later discussions of the accuracy considerations and differential 
methods used in the experiment. 

In the hyperbolic mode, geographic position (GP) is determined from 
time difference (TD) readings provided by a Loran-C receiver. The 
difference in reception time of pulses received from the master station 
and pulses received from each secondary station give the receiver's 
position in terms of a hyperbolic lattice system. Geographic position 
can be calculated from the TD readings using Loran-C lattice tables 
published by the Defense Mapping Agency or computed directly using 
various computer techniques. Nautical charts overprinted with Loran-C 
TD lattices are also available for plotting geographic position 
directly without converting first to latitude and longitude. 

The hyperbolic lattice system is a series of hyperbolas. A pulse 
(actually a group of pulses) is first transmitted by the master 
station. The wave travels to each secondary station of the Loran-C 
chain, and each secondary transmits its own pulse after a time delay 
(called the coding delay) which is unique to each secondary. 
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The master pulse is received first by the receiver. A short time 
later the pulses from the secondaries are received. The time dif- 
ference between reception of the master and each seco n dary is measured 
in the receiver. Each TD identifies one line of position (LDP) in the 
set of LDP's for a master-secondary pair. Figure II-l shows a typical 
Loran-C master-secondary pair with a set of hyperbolic LDP's. Each LDP 
in this set represents a line of equal time difference between 
reception of the master pulse and reception of the pulse from the 
secondary. Figure II -2 shows an example of a hyperbolic lattice system 
formed by two master-secondary pairs. 

As with latitude and longitude, two LDP's are required to define a 
geographic position. Unlike latitude and longitude, two TD readings 
from different secondaries do not uniquely define a GP, since the TD's 
may intersect at two points within the coverage area. This is normally 
not a problem, since even an approximate knowledge of the geographic 
position of the receiver will allow determination of the correct 
intersection point. 

When more than two LDP's are available to define a GP (when signals 
from more than two secondary stations can be received) , the 
intersection of the three LDP's define three geographic points (see 
Figure IV-2) . The accuracy of the two best pairs of LDP's is 
sufficient for most navigation purposes, but it is possible to combine 
three or more LDP's to estimate the position of one point. A method 
using corrections from three LDP's is used in the experiment and will 
be discussed in Section IV. B. 
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FIGURE II-l. Typical set of hyperbolic LDP's associated with a roaster- 
secondary pair. 




FIGURE II-2. Typical hyperbolic lattice formed by sets of LDP's 
associated with two roaster-secondary pairs with a 
common master station. 
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B. Time Difference Equations 



The difference between arrival time of the signal from the master 
and arrival time of the signal from a secondary can be described by 
the following equations (Ref. 1) : 

The basic time difference equation for Loran-C is: 

TO = ts - [II-l] 

where: 

TD = time difference at a given point P. 
tg = time of reception of secondary pulse. 

= time of reception of master pulse. 

A more useful form of this equation takes into account the 
components which make up the travel time of the signals over each 
section of the paths from the stations to point P: 

TO = (T s + DT S + ED) - (T ra + DT m ) + ASF [II-2] 

where: 

TD = Time Difference at a given point P. 

T s = Distance in time units between the secondary station 
and point P. 

= Distance in time units between the master station 
and point P. 

DT S = Secondary Phase Correction for all seawater path 
between secondary station and point P. 

DTjj, = Secondary Phase Correction for all seawater path 
between master station and point P. 

ED = Emission Delay for the secondary station. 

ASF = Additional Secondary Phase Correction for point P. 
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Equation II-2 contains two corrections which are applied to Loran-C 
TD readings. The velocity of electromagnetic waves in air under 
standard conditions is changed when the wave propagates over a surface. 
The Secondary Phase Correction Factor (DT S and ETT m in equation II-2) is 
a time correction for the change in velocity caused by propagation of 
electromagnetic waves over seawater. This change in travel time can be 
calculated by the following empirically derived equations (Ref. 2) : 

For distances greater than 100 statute miles (T > 537 /rsec) : 



where: 

DT = Secondary Phase Correction for all seawater path. 
Aq = 129.04398 Msec 2 

A-l = -0.40758 nsec 

A 2 = 0.00064576438 

T = Travel time (T s or T m in equation II-l) in Msec. 
For distances less than 100 statute miles (T < 537 nsec ) : 



where: 

Bq = 2.7412979 nsec 2 

B;l = -0.011402 /xsec 

B 2 = 0.00032774624 

T = Travel time in Msec. 



DT = Aq/T + + A 2 *T 



[II— 3] 



DT = Bq/T + B;l + B 2 *T 



[II-4] 
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In equation II-2 the Additional Secondary Fhase Factor (ASF) 
corrects for propagation over land (actually any other non-seawater) 
surfaces, which have a wide variety of conductivities (see Section 
III-C-2 or Ref. 3 for further explanation of ASF Corrections) . The ASF 
correction is in addition to the Secondary Fhase Correction for all 
seawater path. 

Secondary Phase Corrections for all seawater path are always added 
in before calculation of the geographic position from TD's. ASF 
corrections cure much smaller, and are often determined from tables 
which require a dead reckoning position within 10' of latitude and 
longitude of the final position. If no DR position is available, the 
position calculated from TD's corrected only for all seawater path is 
used as a DR position to determine ASF corrections . When used in 
hyperbolic mode, the ASF correction in Equation II-2 is a net 
correction which applies to the particular master-secondary pair (the 
difference between the range-range corrections for each station) . 

Emission Delay (ED) in Equation 1 1-2 includes a coding delay, which 
is the time between reception of the master pulse at the secondary 
station and transmission of the secondary's pulse, plus one way travel 
time from the master to the secondary, including the Secondary Fhase 
Correction for all seawater path for the master-secondary baseline path 
(Ref. 4) . 

The result of application of these corrections to the basic time 
difference between reception of master and secondary signals is a 
reduction of the time difference to a path through air with an Index of 
Refraction of 1.000338 (Ref.4:p.4). 
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C. Differential Oorrecticns 



Differential techniques can be used to improve the absolute 
accuracy of Loran-C TD readings. Differential Loran-C uses monitor 
receivers located at fixed points in or near the area where improvement 
is needed. The most common method is to use TD readings of the 
monitor receivers to determine the magnitude of the variations in the 
TD's from the long term average of the TD's at the fixed points. These 
variations are then applied to receiver readings in the area of 
interest before geographic positions are calculated from them. 
Corrections are updated periodically as explained later in this 
section. (Ref.l:p.434) 

The differential correction (including geographic correction) to 
Loran-C readings at an unknown point P 2 with a monitor receiver at a 
known point Pj can be described by the following equation: 

TD at P 2 = Received TD at P 2 + Differential Cor. + Geographic Cor. 

or 

td 2C = td 2 + C^l - TDic) + TDg [II-5] 

where 

TD 2 C = Corrected TD reading at Point P 2 (the unknown point) . 

TD 1C = Long term average TD at P-l (location of the monitor) . 

TDi = TD reading of receiver at the known point (P^) . 

TD 2 = TD reading of the receiver at the unknown point (P 2 ) . 

TDg = TD variation which is a function of location. 
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For practical purposes it is possible to use any correction from 
real time to a correction averaged over a period of 15 ) to several days. 
The method used in this experiment (described in Section IV-C) was to 
determine a regression line of monitor receiver TD readings at the 
known point. A computed correction was applied to the TD reading at 
the unknown point based on the time of the TD reading at the unknown 
point and the equation of the regression line. This method was used 
because only one monitor receiver was available and it was used as both 
a monitor at the known point and at the unknown points. In this 
experiment a correction based on reception of three secondary stations 
(rather than two) was applied to the TD readings from each secondary 
station before the regression line parameters were calculated, as 
explained later in Section IV-B. 

The geographic correction - that part of the difference in TD 
readings at points P! and P 2 which is a due to the differences in ASF 
corrections at the two points - is a particular problem. For 
geographic points at sea, tables of ASF corrections within 10' of 
latitude and longitude cure available, and may be accurate enough to 
provide corrections within the necessary accuracy. Tables of 
corrections for geographic points on land are not published, but were 
obtained for this experiment from the Defense Mapping Agency 
Hydrographic/Topographic Center. 

Test results of differential Loran-C on the East Coast U.S. chain 
have shewn accuracies of about 4.6 m CEP at distances up to 140 km from 
the monitor site (Ref . l:p.434) . 



11 



IH. GENERAL ANALYSIS OF ERRORS 



To quantify the amount of error and compare the accuracy of 
different types of corrections made to Loran-C geographic position 
fixes, it is necessary to provide seme basic definitions. The 
following sections give the definitions used in this study for the 
types of accuracy and measures of accuracy. One and two dimensional 
sources of error in Loran-C positioning are listed and approximate 
values for the error are given where possible. Error reduction based 
on differential correction is described. 

A. Types of Accuracy 

There are several ways to specify the accuracy of electronic 
navigation systems. The following definitions of three types of 
accuracy have been found useful in this study (from Ref. 5): 

Absolute Accuracy : For purposes of this study, absolute accuracy 
is a measure of the ability to determine geographic position (in terms 
of latitude and longitude or other fixed coordinate system) from a 
given set of Loran-C coordinates (Ref .5:p.l72) . The error includes all 
types of systematic and random error associated with the system. 

Absolute accuracy in this experiment refers to the accuracy of the 
calculation of the positions of the "unknown" points based on receiver 
TD readings, three station corrections, and differential corrections, 
compared to the known published positions of those points. 
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Repea tabl e Accuracy : The carmen navigation definition of 

repeatable accuracy is the ability to return to the same point using a 
given set of electronic coordinates (Ref.5:p.l72) . Error is measured 
in linear units. This differs from absolute accuracy in that the 
position itself need not be kncwn in terms of scare other coordinate 
system. Repeatable accuracy is normally better than absolute accuracy 
because the part of systematic error which is a function of position 
has been removed (the position remains the same) , leaving only random 
error and time dependant systematic error. Repeatable accuracy in the 
context of this experiment is a measure of the precision of measurement 
of the TD's, which is most affected by the electronic circuitry of the 
transmitters and receivers involved (over short periods of time in 
which terrain and atmospheric factors do not change significantly) . 

Repeatable accuracy is similar to the accuracy expected from 
differential corrections. The difference between the two is that 
since repeatable accuracy refers to the ability to return to the same 
point using the same set of received TD readings, there is no 
difference in the ASF (terrain) correction but there may be some 
difference due to variability of TD's with time. 

Differential corrections have different ASF corrections 
at different points. These must be calculated and applied 
to the TD readings. At a single point the ASF corrections 
will probably not change over the time period of interest, 
but over longer periods some error will be introduced. The 
uncertainty in ASF corrections could be a significant 
source of error. 
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Relational Accuracy : The accuracy with which it is possible to 
determine the position of one observer in relation to another observer 
using the same system (Ref .5:p.l72) . A form of relational accuracy 
enters into differential corrections when two different receivers are 
used. Even at the same point at the same time there will usually be a 
slight difference in the loran-C readings of different receivers. This 
is caused by errors associated with the receiver circuitry and random 
errors over the slight differences in the time at which the readings 
are made by each receiver, since all other factors are the same. The 
same type of error exists with two different receivers separated by 
some distance (as with differential corrections) , although this type of 
error is usually small compared to other sources of error. 

B. M easures of Accuracy 

In addition to the types of accuracy listed above, it is necessary 
to specify how accuracy is to be measured. This can be viewed as an 
estimate of the amount of error. The following are definitions of the 
basic methods used here of specifying the measures of accuracy of 
position. Their usefulness and limitations in this application are 
also described. 

Root Mean Square (RMS) Error : This is the linear (one dimensional) 
error in the location of an LDP (Ref .5:p.l72) . It can be specified in 
linear units or time units (meters and microseconds are used here) . 
The normal specification of RMS error is the number of meters or 
microseconds representing sane whole number of standard deviations from 
the mean. In the application used here, the variation in TD readings 
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fron each master-secondary pair can be expressed as a standard 
deviation in microseconds, which corresponds to a displacement of the 
LDP by a corresponding number of meters measured normal to that line of 
position. 

Circular Error Probable (CEP) : This is the two dimensional error 

which specifies the radius of a circle within which there is a 50% 
probability of being located, based on positions calculated from the TD 
readings of the receiver. The estimated position of the point (based 
on latitude and longitude or other coordinate systems) and the radius 
of the circle are used to specify this measure of accuracy. Circles 
other than the standard 50% circle can be used (67% and 90% are 
common). (Ref.5:p.l73) The main advantage of using CEP (rather than 
the error ellipse described below) is that direct comparison between 
accuracies of different systems, techniques, or at different points can 
be made. CEP in this experiment is used as a measure of the precision 
of the TD's indicated by the receiver at a fixed point and the 
precision of corrected TD's based on these uncorrected TD's. 
Corrections for the absolute position (the bias of the TD's from true 
readings) is of secondary importance in this application. 

Circular error estimates in this experiment are based 
on the radius which 50% of the data are greater and 50% are 
less, as calculated by a sorting program. 
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Error Ellipses : Error ellipses are the preferred way of 



describing the accuracy of electronic navigation systems for many 
purposes. As with circular error, accuracy is specified in terms of 
the probability of an estimate of geographic position falling within a 
geometric figure centered on the actual geographic position (an ellipse 
rather than a circle is used) (Ref. 6) . Unlike the circle used in CEP, 
an ellipse provides additional information about the directional 
characteristics of the accuracy, which can be very important in some 
applications. The specifications for elliptical accuracy are usually 
the lengths of the major and minor axes and the orientation of the 
ellipse (such as the azimuth of the direction of the major axis) . A 
disadvantage of using error ellipses in this experiment is the 
difficulty of directly comparing the accuracy of different ellipses 
with different orientations. For that reason most of the accuracy 
comparisons in this study have been specified in terms of circular 
error. Error ellipses have been used where appropriate . 

Error ellipses are related to the orientation of two intersecting 
LDP's by the following equations (Ref .6:p. 72-74) : 



a 2 2 sin 20 



tan 2© 



[III-l] 



a 1 2 (cos 20) + a 2 2 





[III-2] 
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[HI - 3] 




°1 2 + °2 2 ~ ( ( a l 2 + a 2 2 )^ ” 4 (sin20)a 1 2 a 2 2 ) ^ 



2 sin 2 <p 



where: 

0 = the angle of the major axis of the ellipse 

from the direction of the LOP with smaller 
variance, toward the LOP with greater variance, 
in the direction of the smaller angle. 

Oi 2 = variance of the LOP with smaller variance. 

a 2 2 = variance of the LOP with the greater variance. 

0 = the smaller angle between LOP's. 

These equations were used to calculate the dimensions and 
orientation of the error ellipse for Range-7, the "known" point used 
for differential corrections (see Figure V-3) . 



C. Sources of Linear Error in Loran-C Time Differences 

The position fixing accuracy of the LDRAN-C system depends on a 
number of factors. These can be divided into two categories 
(Ref .7:p. 18) : temporal errors affecting the repeatability of the 

system (errors which vary with time) , and systematic (time invariant) 
errors. Both of these types of error nay also be functions of 

'geographic position. 

The effect of these errors is to produce TD readings at the 
receiver which differ from those which would be expected on a perfect 
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hyperbolic lattice on the surface of a perfect ellipsoid, resulting in 
error in the calculation of the geographic position. A linear error of 
1.0 /nsec in time difference is approximately equivalent to 150 m on the 
baseline, with one dimensional error at other points being 15 ) to about 
ten times greater, depending on the location within the geometry of the 
system. 

The main sources of error in the determination of TD's used to 
calculate geographic position can be classified as follows (Ref. 8 ) : 

• Transmitter synchronization timing error's. 

• Errors associated with propagation effects. 

• Atmospheric noise. 

• Receiver errors . 

1. Transmitter Synchronization Timing Errors 

Since time differences between reception of a master pulse and 
secondary pulses are used to calculate geographic position in 
hyperbolic Loran-C, it is essential that the synchronization of timing 
between transmissions from the master and secondaries be highly 
precise. 

In order to accomplish this, each station is individually timed 
using cesium clocks which drive pulse and group repetition interval 
timing circuits. Independent monitor receivers are used to provide 
data to computers which adjust the relative timing between master and 
secondaries of the chain in small steps (usually 10 to 20 ns) . 
Absolute timing of the master stations is adjusted to operate in 
coordination with Universal Time, Coordinated. (Ref. 8 :p. 1129) 
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Estimates of timing error due to transmitter synchronization 
are on the order of 0.05 jisec (Ref .9:p. 537) . Actual data from an 
operational Loran-C chain shewed time differences controlled within 
0.024 /isec of the standard value, with a standard deviation (in 
absolute time) of 0.032 /isec (Ref. 10: p.227) . A transmitter 
synchronization error of 0.01 /isec is believed to be attainable if 
ground monitor data is used to remove the mean synchronization error 
(Ref. ll:p.401) . 

2. Errors Associated with Propagation Effects 

This includes most types of error external to the electronics 
of the system, and is the main source of error in the determination of 
geographic position. 

Factors which influence groundwave propagation include the 
following (Ref .12: pp. 39-53 and Ref. 13 :pp. 173-187) : 

• Electrical Conductivity of the surface material 

• Dielectric Constant of the surface 

• Index of Refraction of the air 

• Lapse Rate of the Index of Refraction with altitude 

These factors are complex functions of several variables, the 
main ones being: 

• type and composition of the surface over which the wave travels 

• topography of the surface 

• meteorological conditions 
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At any particular geographic position the effect of these 
variables may or may not be considered a function of time, depending on 
the time periods of interest. 

Also of importance is the effect of these factors on TD as a 
function of geographic position. For each different geographic 
position, the electromagnetic waves from each transmitter must traverse 
a different route, with different types of surface, topography, and 
atmospheric conditions, to reach the receiver. 

The effects of irregular terrain are practically time 
invariant, while ground conductivity and dielectric constant may vary 
seasonally (due to changes in moisture content of the soil, changes in 
vegetation, or other factors) . Observed cyclic annual variations on 
the order of 0.5 ju sec are probably due to seasonal variations in ground 
conductivity and dielectric constant and to seasonally varying 
atmospheric conditions. 

Meteorological conditions sufficient to cause noticeable 
variations in TD readings may change in very short time periods (days 
or hours) . Short term variations due to meteorological effects (cold 
front and warm front passage through the line of propagation) on the 
order of 0.025 nsec have been observed. These correlate with changes 
in the dry adiabatic lapse rate of the index of refraction . 
(Ref. 14 :p. 1119) . 

For most differential Loran-C purposes it is necessary to 
consider only those temporal propagation effects which occur in short 
time periods (hours or days at most, for most applications) . The 
variation of propagation effects with geographic position must also be 
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considered for precise positioning. Many of these variations cam be 
corrected by a computation of the Additional Secondary Phase Factor 
(ASF). 

The Additional Secondary Phase Factor (ASF) is a measure of the 
systematic error resulting from differences in velocity of electromag- 
netic waves over ground with different conductivities and dielectric 
constants. The Defense Mapping agency vises a computerized system of 
maps with estimates of ground conductivity (effectively including 
dielectric constant) over large areas, and calculates ASF corrections 
in the Loran-C coverage areas. These corrections are published as 
tables of ASF corrections to be used in correcting Loran-C TD readings 
from receivers. 

ASF corrections can be as large as 4 microseconds (Ref. 3). 
The corrections provided in the Defense Mapping Agency tables require 
estimates of the geographic location of the receiver to the nearest 10' 
of latitude and longitude. ASF corrections based on the kncwn 
positions of the survey points used in this experiment were provided by 
EMA using a computer program which calculates the corrections using the 
Millington Method (Ref.3:p.V). 

3 . Atmospheric Noise 

The effects of atmospheric electrical noise are reduction of 
the range at which Loran-C signals can be effectively received and an 
increase the error in TD's computed by receivers in marginal reception 
areas. This is primarily dependent on the circuitry of the receiver. 
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In marginal areas the effect of atmospheric noise will be an 
increase in the standard deviation of the TD's and increased difficulty 
in locking on and maintaining reception of the weaker signals. 
Atmospheric noise varies with time and geographic location, with a 
significant diurnal variation (noise levels are higher at night) . In 
this experiment, significant effects due to local sources of electrical 
noise were noticed in the Monterey harbor area. 

4 . Receiver Errors 

Receiver errors are of two main types: 

• timing error caused by 1 imitations of the internal clock 

circuitry of the receiver 

• errors caused by processing time of the receiver 

For a receiver at a fixed geographic location, processing time 
is not a source of error. The receiver will be at the same location 
regardless of the time required to output the TD's or geographic 
coordinates. The main source of receiver error at a fixed location 
(other than the noise considerations discussed above) is the timing 
error in the internal clock circuitry. Loran-C receivers vary greatly 
in their timing errors. The Racal Megapulse Monitor receiver used in 
this experiment has a specified timing error of 0.02 fisec or better. 

For mobile receivers the processing time of the receiver 
circuitry can be important. After the signal is received it must be 
processed and either stored in the form of electronic information or 
visually displayed. Error is introduced if the receiver moves during 
this time interval . 
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Obviously, the amount of error in the determination of 
geographic position from this source depends on the processing time and 
the speed of the receiver. Receivers which output data in the form of 
geographic coordinates must calculate latitude and longitude from the 
TD readings of the receiver, requiring more processing time and 
resulting in greater error fron this source. 

At the speeds normally used in hydrographic surveying 
(typically less than 20 kt) the processing time is insignificant for 
practically all receivers. For data which is stored in the form of 
receiver output TD's for later correction and conversion, as with the 
Racal Megapulse receiver used here, error should be insignificant even 
at much higher speeds. 

D. Geometric Dilution of Precision 

This source of error in geographic position results from the 
crossing angle of the Loran-C LDP's (lines of position). Figure IV- 1 
shews two hyperbolic LDP's with standard deviations and a 2 (in time 
difference) and crossing angle a. The optimum crossing angle of any 
two LDP's is 90°, resulting in minimum error in geographic position for 
a given set of standard deviations. 
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E. Summation of Errors 



1. Absolute Accuracy vs. Precision 

The absolute accuracy of an electronic positioning system can 
be viewed as a measure of its ability to give the geographic position 
of the receiver. Precision is a measure of the variation of individual 
time difference readings. If the receiver is set up at a point for a 
sufficient amount of time, many of the random errors can be averaged 
out; the average of a large number of readings can provide a basis for 
determining absolute accuracy. 

In practice this is more difficult than it would seem, since 
seme factors which do not vary significantly over short periods of time 
(hours or days) may vary significantly over longer periods (weeks or 
months) . It is necessary to consider the time periods involved in 
order to determine which factors to consider random and which to 
consider "bias" in the system. 

For this experiment, the sources of random error are considered 
to be transmitter timing errors, seme meteorological conditions, 
atmospheric noise, and most receiver errors. Errors caused by long 
term "bias" of the system are due primarily to ground conductivity and 
type of surface, seme meteorological effects, and any non-random bias 
in the timing of the transmitters. 

Taking these factors into consideration, data taken at the 
"known" point (Range-7) of the differential system was processed to 
remove random error by averaging over a period of time, while long term 
error was removed by comparison of the TD's to their expected values 
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derived from calculations provided by computer computations of the TO 
at that point. 

Errors which may change within the time frame of the experiment 
were reduced by a linear regression of the trend of data taken at the 
known point and projected to the correction at the unknown point 
mathematically . 
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IV. THREE STATION AND DIFFERENTIAL CORRECTIONS 



A. Basic Geometry of Ioran-C Hyperbolas 

Ibis section provides a summary of the information necessary for 
the analysis of Loran-C corrections for three secondary stations and 
for differential Loran-C corrections presented in parts TV.B and IV. C, 
adapted from Laurilla (Ref.l) and others as noted. 

1. Lanewidth 

The baseline of a Loran-C master-secondary pair is the shortest 
path over which the transmitted signals may travel between the master 
station and the secondary station. On a spherical earth model it is 
the arc of the great circle through the positions of the master and 
secondary stations. On an ellipsoid it is the geodesic from the 
master to the secondary station. In the set of hyperbolas of time 
difference for the master-secondary pair, the separation between 
hyperbolas is a minimum constant value along the baseline. 

A lane is the distance between two adjacent hyperbolas of unit 
time difference. The unit normally used is the jizec (10“ 6 sec). In a 
series of hyperbolas, as shewn in Figure II-l, it is obvious that the 
lanewidth varies with location. At any point, lanewidth can be 
viewed as the inverse of the gradient of the ratio of the change in 
time difference to the change in distance (Ref. 15), or: 

L = { | (dTD/dl) | r 1 [IV-1] 
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Figure IV-1. Typical intersection of two lines of 
position. Most probable point is at the intersection of 
the two LOP's. 




Figure IV-2. Typical intersection of three lines of 

position. When three or more secondary stations can be 
received at a geographic position, calculated position 
depends on the LOP's which are used to make the 

calculation. 
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The lanewidth at any point in the set of hyperbolas can be 
expressed in terms of the lane width on the baseline. Baseline lane 
width can be expressed by the following equation: 

^ = -f“ [IV-2] 

where: 

1$-, = Lanewidth on the baseline (meters//xsec) 

c = Speed of electromagnetic waves on the baseline (m/^isec) 
(note that time units are in the unit lanewidth) 

Lanewidth at any other point (not on the baseline) can be 
expressed by equation IV-3a, where /3 is the angle subtended by the 
baseline from the point. Note that this angle is the absolute 
difference between the azimuths from the point to the master and from 
the point to the secondary station. 

L = I b ( sin 03/2) ) -1 [IV-3a] 

or 

L = IfeG 

where: 

L = lanewidth at a point (per unit of time) 

I b = Lanewidth on the baseline (per unit of time) 
f3 = Angle subtended by the baseline from the point 
G = The Lane Expansion Factor, G = (sin 03/2) )“■*- 

Note that the lane expansion factor depends only on the angular 
separation of the master and the secondary stations. This obviously 
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depends on the distance of the point frcsn the baseline center and the 
angular separation of the baseline and the vector frcsn the baseline 
center to the point. 

Equation XV-3a provides a means of calculating the normal 
distance between any two hyperbolas, provided they are near enough that 
their lanewidths are not significantly different: 

d = Lt [IV-3b] 



where: 

d = distance in meters between the two hyperbolas. 

L = lanewidth at the point of interest ( in m/iisec ) . 

t = time difference between the two hyperbolas in jxsec. 

2. Direction of Hyperbola 

The direction of a hyperbola at any point on the hyperbola is 
along the tangent at that point, directed away frcsn the baseline of the 
master and secondary stations used to define the hyperbola. For 
applications in this experiment, Laurilla (Ref 1) provides the follow- 
ing useful definition of the direction of a hyperbola: "... the 

direction of a hyperbola at any point coincides with the bisector of 
the angle formed by the lines joining the point to the pair of 
stations." (Ref.l:p.92) 

For purposes of this experiment, it is convenient to express 
the direction of the hyperbola by its azimuth. This azimuth can be 
expressed in terms of the azimuths (from north) of the master and 
secondary stations by: 



29 



“h 



(ctm + a s )/2 (where 1 0 ^- 05 ! >180*) 



[IV-4a] 



oth = 180* + (<% + Q^)/2 (where lom-Os^lSO”) [IV-4b] 

where: 

= azirmith of the direction of the hyperbola, 
ctjn and a s = azimuths of the master and secondary from the point. 



3. Applications of Geometry of Hyperbolas to Error Analysis 

The geometric description of loran-C hyperbolas in parts IV. A. 1 
and IV. A. 2 is particularly important in the analysis of error and the 
derivation of corrections used in this experiment. Some of the direct 
applications are outlined in the following paragraphs. 

The error expressed as an error ellipse requires the directions 
of hyperbolas in order to determine the orientation of the ellipse. 
The dimensions of the ellipse require the conversion of variances of 
time differences of each hyperbola to variances of distance in the 
major and minor axis of the ellipse. Standard deviations of time 
differences can be converted to distances by the following equations: 

cr d = cr-t-L [IV- 5a] 
a d 2 = a t 2 L 2 [IV-5b] 



30 



The correction for three secondary stations uses variances in 
time differences converted to distances as an input, using equation 
IV-5b. Converting time difference to distance in this way correctly 
weights the precision of the LDP's when determining the correction 
factor based on geometric considerations. 

B. Correction for Three Secondary Stations 

In situations where three secondary stations can be received, the 
result is three LDP's which do not, in general, intersect at one point 
(see Figure IV-2) when they are used to calculate the geographic 
position, although they must in fact intersect at the common point at 
which the readings were taken concurrently. Under usual navigation 
conditions the two best LDP's are used to determine a geographic fix. 
A drawback of this method is that it fails to consider all of the 
available information (i.e. one of the available lines of position is 
not used in the calculation of the geographic position) . By utilizing 
all three LDP's it should be possible to improve the accuracy of the 
geographic fix. 

Sections B.l and B.2 will develop a mathematical adjustment which 
can be used to calculate the most probable corrected TD readings of a 
Loran-C receiver used as the fixed receiver in a differential Loran-C 
system utilizing three LOP's (i.e. receiving the master and three 
secondary stations) . When the correction is applied, the three 
corrected TD's will intersect at the most probable point, based on long 
term average values of the TD's. The derivation of this adjustment 
could take any of several possible forms, but a primary cons iderat ion 
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here is the need for a form which is easily applicable to computer 
processing. Section B.3 describes a method of applying the corrections 
for non-differential use, using only one receiver and with no knewn 
fixed pxoint. 

Section B.4 will develop regression equations for calculating the 
corrected TD's directly from time for use as a differential correction 
for a remote receiver. Note that the methods developed in these 
sections are primarily for use in the differential Loran-C computations 
necessary for this experiment. 

1. Geometric Considerations Used to Develop a Correction Factor 

The following assumptions are made for this derivation: 

• Each LDP has a unique direction (described in Section IV. A. 2) 

• Long term average TD's (designated TDl av , TD2 av , and TD3 av ) for 

for the known point are available 

• Variances of the TD values (designated a^ 2 , a 2 2 , and a 3 2 ) are 

known or can be calculated 

• Variations of TD's from their average values is in the form of a 

Gaussian distribution 

• Over the relatively small area being considered, the surface can 

be considered a plane and LDP's for the same master-secondary 

pair can be considered parallel straight lines 

Several things should be noted about these assumptions. Long 
term averaging times for this experiment will be on the order of a few 
hours during acquisition of data; for purposes of the derivation, 
standard deviations are converted between time units (/xsec) and 
distance units (m) using Equation IV-5a; and, the assumptions of a 
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N 




FIGURE IV-3. Geometry used in the derivation of a correction factor 
for three LDP's not intersecting at the same point. 
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Gaussian distribution, plane surface, and parallel LDP's are reasonable 
for this type of data (Ref. 6 :p. 8 ). When two LDP's with different 

standard deviations intersect (see Figure IV-1) , the point of 
intersection is the most probable geographic position (Ref . 6:p.8-12) . 
The probability of the location of the true position is best described 
using an error ellipse as described in Section III-B. 

A typical situation is shown in Figure 3, a set of three 

concurrent LDP's intersecting at three different points. Point P is 
the location of the monitor receiver used to determine differential 
corrections, with TDl av , TD2 av , and TD3 av being the long term average 
TD values at that known point (note that in figure TV-3, the LDP's of 
TDl av , TD2 av , and TD3 av are not shown, but are parallel to TD^, TD 2 , 
and TD 3 , respectively, and pass through point P, as shown in Figure 
IV-4 ) . Although the long term average LDP values intersect at point 
P, triangle Pi2 p i3 p 23 represents a triangle of fixes formed by the set 
of TD values which may occur from any concurrent set of TD's, which are 
unlikely to equal the long term average. Each set of TD readings 
(which occurs in the data every few seconds) will result in a new 
triangle of fixes. 

The vectors r^, r 3 , and r 3 represent the difference between the 
individual TD values in a concurrent set of data and their respective 
long term TD values. The r vectors are at right angles to both the 
long term and current LDP's. For example, if the long term average of 
TD-l (TDl av ) is 29939.200 and the value of TD^ in a set of TD readings 
currently being considered is 29939.300, the value of r 3 for the 
current set of TD readings is 0.100 (all in ^isec) . When converted to 



34 



meters using equation IV- 3b, this is the distance between the average 
TD values (LOPi av ) and current TD values (TD-^) . These LDP's can be 
considered as parallel straight lines over the relatively small 
distances considered for this correction. Note that there are three 
LDP's in each set of concurrent data representing TD values for three 
master-secondary pairs. The master station is common to each of these 
pairs. 

Within the triangle of fixes Pi2 p i3 p 23 fanned by the three 
LDP's of a set of data, point P' represents the most probable point of 
the true geographic position indicated by that particular concurrent 
set of data compared to the long term average data. The distances of 
P' from each LDP (the lengths rp' ) are related to the relative 
precision of each LDP, P' being closer to LDP's with higher precision 
(to be discussed in more detail in the next section) . 

The problem must be solved in general terms so that any set of 
TD's in the data can be corrected to a new set of TD's which intersect 
at the most probable point. Since the triangles formed by all sets of 
data are similar triangles, the lengths of q 1 , r 2' ' anc ^ r 3 ' > will have 
the same proportion for each set. The corresponding r^' values between 
sets of data will be related by a common factor, which will be 
designated as the Correction Factor, k. 

Other information which should be noted about Figure XV-3 
include the follcwing: 

• The following sets of lines can be considered parallel over the 
area represented by the figure: TDl av || TD^ ; TD2 av || TD 2 ; 

TO3 av II ; r l II r l ' ; r 2 f r 2 ' ; and r 3 || r 3 ' . 

(continued on next page) 
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• The following sets of lines are perpendicular: and rj_' J. 

TDl av and TD^ ; r 2 and r 2 ' ± TD2 av and TT^ ; r 3 and r 3 ' J- 
TD3 av and TE^. 

• The azimuths (from north) of the directions of the LDP's 

represented by TDl av , TD2 av , and TD3 av are og, a 2 , and a 3 , 
respectively. These are the same as the azimuths of TD 1( Tt^, 
and TL^, respectively. 

• The following angles of triangle Pi2 p 13 p 23 are defined by: 

a 12 = a 2 “ a l » a l3 = ff/2 ~ (a 3 - «i) ; and a 23 = a 3 - a^. 
These angles will be the same in the similar triangles formed 
by any lengths of r^, r 2 , and r 3 . 



• LDP's 1 and 3 have the largest absolute angle (+ or -) between 
them, and LDP 2 is the intermediate angle. 



The next section will provide a mathematical derivation of the 
expression for a correction factor to determine the lengths of r^ ' , 
r 2 ' , and r 3 ' and the corresponding time adjustments to TD^, TD^), and 
TI^, based on the geometry presented in this section. 



2. Mathematical Derivation of the Correction Factor 



An important question in this analysis is how to weight the 
relative precision of each LDP in such a way that it can be used in 
the calculation of a correction factor. Reference 16 shews that for 
this type of data the "weights are inversely proportional to variance" 
(Ref.l6:p.60) , meaning that the lengths of the r^' corrections should 
be proportional to the variances (in units of length) of the 
corresponding TD values for each LDP. Using units of length, in the 
general solution for all sets of data, the following equation applies: 

rj/Li = ka^Lj 2 [IV-6] 

where : 

r^' = corrections for TD's 1, 2, and 3 (in /i sec). 
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o± 2 = variances of TD's 1, 2, and 3 (in /isec 2 ) . 
k = the Correction Factor (units are meters) 

= lanewidths of TD's 1, 2, and 3 (in ia/fis ec) . 

Note that a^ 2 L^ 2 is a unitless weighting factor. 

For any triangle there will be seme unique value of k which 
relates the distances from point p' to each LOP to the variances of the 
LOP's according to equations IV-6. To obtain a general expression for 
k which can be easily calculated by computer, it will be helpful to use 
only the known values of the TD's, variances, azimuths, and lanewidths 
for each LDP. The r values are readily obtained from each set of TD 
data by subtracting the long term average values from current TD 
values. Values of a^ 2 and are constants for this small area. 

If each LOP is considered a straight line and a cartesian 
coordinate system is assumed with north as the +y and east as the +x 
direction of the coordinate axes, then each LOP can be specified by the 
general equation for a straight line: 

y = mpe + b^ (i = 1, 2, or 3) [IV-7] 

Three such equations can be written (one for each LOP) , giving 
three linear equations in two unknowns, x and y. These equations will 
not have a common solution when the three LOP's do not intersect at a 
common point. 

For each master-secondary pair, all LOP's of the set will have 
the same slope (m^) within the small area being considered (TD^ is 
parallel to TDi av for any TDjJ . The y-intercepts (the b^'s) will be 
zero for the average LOP's intersecting at point P. Y-intercepts for 
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other LDP's with non-zero r values can be calculated from the 



relationships between azimuths (denoted ajJ and the r values of each 
LDP by the following equation: 

IT ' 1C * 

bi = = — (i = 1, 2, or 3) [IV-8] 

cos(a-L-7r/2) sin a j_ 



N (y) 




FIGURE IV- 4 . Typical average LDP through the origin of a cartesian 

coordinate system and corresponding parallel LDP of an individual 
TD reading. Slopes are the same for both. Y intercept is the b 
value in equation for the LDP (zero for the average LDP) . 
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Slopes of the TD's can be expressed in terms of cp as: 



m^ = -tan <p± [IV-8b] 

Figure IV-4 shews the assumed coordinate system along with the 
LDP's for which the slopes and intercepts are described. When 
corrections are applied, the corrected r values are as follows: 

r iC L i = r i L i “ to i 2 Li 2 [IV-9] 

Note that in the equation for corrections (IV-9) the units of 
r^oL^ are meters. If a^Lj 2 is considered a unitless weighting factor, 
then the units of k must be meters. This agrees with the units of 
Equation 13. 

The value of r^ is calculated by: 

*1 = TDi - TDi av [IV-10] 



where: 

r-[ = deviation of an LDP from the long term average (nsec ) . 

TE>i = time difference reading of a master-secondary pair at 
point P. 

TDi av = long term average of the TD for the same pair at 
point P. 

Figure IV- 3 shews that the r^' corrections must be subtracted 
from r^ and r 3 , and added to r 2 , for the configuration shown. It can 
be demonstrated that the r^' corrections are always of the same sign 
(both + or both -) for the two most widely spaced LDP's and of the 
opposite sign for the intermediate LDP, where the direction vectors at 
point P are in the positive directions of increase of TD and the angle 
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is in absolute terms (maximum being 180°). This is true for all 
positive and negative r values. Positive and negative k values will 
reverse the signs of all corrections in equations IV-9 so that they 
will be correct for all configurations. 

The corrected r values of equations IV-9 can be substituted 
into equations IV- 8 , providing a set of y- intercepts (the b^ values in 
the set of general equations TV-7) for the corrected LDP's which 
intersect at point P' . These intercepts contain the unknown Correction 
Factor k. Note that the slopes of the corrected LDP's are equal to 
those of the corresponding uncorrected LDP's, meaning that the m^ 
values in equations TV-7 are unchanged. 

Substituting the corrected y- intercept (b^) values into 
equations IV-7 results in three linear equations in three unknowns: x, 
y, and k: 



r^Li - ka^L ] 2 

y = (-tan 4>i)x + [IV-lla] 

cos 02 



r 2 L2 + ka 2 2 Lr) 2 

y = (-tan 0 2 )x + [iv-llb] 

cos 0 2 



r 3 L 3 “ to 3 2h 3 2 

y = (-tan 0 3 )x + — [IV-llc] 

cos 0 3 
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Equations IV-11 can be rewritten as the following set of linear 



equation in the unknowns x, y, and k: 



° 1 2 l 1 2 r lLl 

(tan 0i)x + y + ( )k = 

cos <pi cos <p^ 

° 2 2 l 2 2 r 2^2 

(tan 0 2 ) x + y - ( ) k = 

cos 0 2 cos 0 2 



o 3 2 L3 2 r 3 L3 

(tan 0 3 ) x + y + ( ) k = 

cos 0 3 cos 0 3 



[IV-12a] 



[XV-12b] 



[IV-12C] 



Solving equations IV-12 by Gaussian elimination gives the following 
general expression for the Correction Factor k: 

tan 0 3 - tan 0 2 tan 0 3 - tan tan 0 2 - tan 0^ 

L l r l “ l2 r 2 + i<3 r 3 

cos 0± cos 02 cos 0 3 



tan 0 3 - tan 0 2 tan 0 3 - tan 0 3 tan 0 2 - tan 0^ 

L^a^ 2 + L^ 2 ^ 2 + L 3 2<y 3 2 

cos 0i cos 0 2 cos 0 3 

[IV-13] 

where: 

r^ = differences from average TD's (/isec) . 

k = Correction Factor for three master secondary pairs (m) . 

0j_ = directions of increasing TD (azimuth from north) . 

a^ 2 = variances of LDP's (part of unitless weighting factor) . 

Li = lanewidths (in meters/^sec) 

Note that Li 2 ^ 2 terms in denominator are unitless weighting factors. 
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Equation IV -13 may be simplified by the following substitutions: 



A 1 = L]_(tan 03 “ tan <£2) /(cos 03) [IV- 14 a] 

A 2 = L2 (tan 03 - tan 03)/(cos 03) [IV- 14 b] 

A 3 = L3 (tan 02 - tan 03)7(005 03) [IV- 14 c] 

A4 = L3C7 + L2Cr 2 2 A 2 + L 3 Q 3 ^A 3 [IV— 14 d] 

Resulting in the following equation for k: 

k = (A^t-l - A 2 r 2 + A 3 r 3 )/A 4 [IV- 15 ] 

The x and y coordinates of point P' can be calculated by: 

r^cos 02 - r 2 cos 03 - k03 2 cos 03 - ka 2 2 cos 02 

x = [IV- 16 ] 

(cos 03) (cos 02) (tan 03 - tan 03) 

y = - (tan 03) x + (r3 - a3 2 k)/(cos 03) [TV- 17 ] 

Two other useful parameters can be calculated from the values 
of x and y. These are the azimuth and range of the of point P' from 
point P. These are given by the following: 

Azimuth (from north) : 

a p pi ~ n /2 - tan -1 (y/x) (x > 0 ) [IV- 18 a] 

apPi = 37 T /2 - tan“l(y/x) (x < 0 ) [IV- 18 b] 

Note: when x =0 the azimuth is zero for positive y and n for negative y 

Range (P to P'): R pp i = (x 2 + y 2 ) 5 * [IV- 19 ] 
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Using the computed value of k for a set of concurrent TD data, 
a correction can be made to the TD values using equations IV-9. Since 
IV-9 is in meters it must be converted to nsec using equation IV-3b. 
The resulting corrected TD's will intersect at the point P' . 

Equations IV-15, IV-16, and IV- 17 do meet the requirements 

stated at the beginning of this section. The inputs are limited to 
azimuths, variances, long term average TD values, and concurrent TD 
data, and the values of k, x, y, and the corrected TD values can be 
computed directly from these by equations which can be calculated 
easily by a computer. 

3. Non-Differential correction for Three Secmdary Stations 

When using the method of Section IV. B. 2 with a single Loran-C 
receiver (i.e. not using differential methods) a correction can be 
calculated based on existing computer programs. While not used in this 
experiment, the method of correction will be outlined in this section. 

Since variation from a known point cannot be used, the two best 
LDP's can be used as the basis for this calculation. While not as good 
as using a kncwn fixed survey point, this will provide a correction 
based on three LDP's. The method is outlined below: 

1. Using the two LDP's with the best fix accuracy (available from 

existing computer programs) , calculate the geographic 
position of the point of intersection. This corresponds to 
point P. 

2. From the geographic position of P, calculate the TD of the 

third master-secondary pair (the one not used in the previous 
step) . 

(continued on next page) 
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3 . 



Using point P calculate the azimuths (from north) of the LDP's. 
Number the TD's as described in the previous sections, with 
the intermediate TD as TE >2 and the TD's with the greatest 
angular separation as TD^ and TD 3 . 

4. Use TD variances calculated over short time periods in the 

general area, or use the values of meters/microsecond cal- 
culated using equation IV- 3 as the approximate weights of 
each LDP. Substitute these for the a 2 values in equation 
IV-13 . 

5. Assume r values of zero for the LDP's with the best fix 

accuracy used in step one. 

6 . Calculate the difference between the observed third TD and the 

calculated value obtained in step 2. Use this for the third 
r value. 

7. Calculate corrected r values as in the previous sections. 

8 . Using any two corrected TD's, calculate geographic position. 

This is the position corrected for three secondary stations. 



Note that all of the above calculations are based on computer 
methods and programs already in use (Ref. 2) allowing corrected posit- 
ions to be calculated readily using the method outlined. 



C. Differential Corrections 



One of the main objectives of this experiment is to apply 
differential corrections obtained at a known point to time difference 
data taken at unknown points . One problem encountered was the 
availability of only one Loran-C monitor receiver, making simultaneous 
data collection at two points impossible. Differential Loran-C does 
not require continuous updating of the differential correction for 
general navigation purposes (hourly or daily corrections are possible) , 
but the need for the most accurate corrections possible for the 
experiment required a continuous update of the correction. 
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1. Mathematical Description of the Differential Correction 



Data from the experiment consisted of a series of readings at 
each point. Data was recorded at intervals of 5.0 seconds for much of 
the experiment, although other time intervals were also used. Data 
consisted of time versus TD for up to four possible master-secondary 
pairs and the signal to noise ratios for each, in the following format: 

Time (GMT) GRI TD X TD 2 TD 3 TD 4 S/N x S/N 2 S/N 3 S/N 4 

Only three master-secondary pairs were receivable in the area 
of the experiment. One survey point (Range 7) was designated the 
"known" point, from which differential corrections were calculated for 
data taken at the other points. Over a period of several days, data 
was collected at the known point for periods of 20 to 30 minutes at 
irregular intervals spaced at about one to three hours (data was 
collected at the other points in the intervening periods as described 
in Section A. 2) . 

To provide estimates of the continuous TD readings at the fixed 
point during all time periods (including the periods between data 
collection at that peint) , a least squares curve fitting process is 
used. The result is three equations: TD as a function of time, for 
each master-secondary pair. 

These three equations can be used to provide an estimate of the 
TD values at the fixed px>int at any time. The difference between the 
calculated TD at a certain time and the long term average TD is the 
Differential Correction at that time. The ASF corrections at the fixed 
point and an estimate of the ASF correction at the unknown point (as 
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described in Section II-B) also enter into the Differential Correction 
at the unknown point. An equation is written for each LDP: 

Cd = TD c - TD av + (ASF U - ASF f ) [IV-20] 



where: 

= Differential Correction at the unknown point. 

TD C = TD (at kncwn point) calculated from regression equation. 
TD av = long term average TD at the known point. 

ASF = ASF correction at the known and unknown points. 

2. least Squares Regression Equations for Time Difference 

Linear regression equations for each TD can be written in the 
following form (Ref. 17): 



TD = b 0 + b x t [IV-21] 

The b-[ coefficients of the regression equation can be found by 
solving the following set of linear equations: 



nb 0 + b]Et^ = STD^ 

b()2t-[ + bjEtj 2 = EtjTDi 



[XV-22a] 



An alternate method of calculating the b^ values directly is: 



nStiTDi - (Sti) (ZTDj_) 

b]_ = 

nStj 2 - (Eti) 2 

ETD^ ETD^ 

b 0 = - b 2 

n n 



[IV-22b] 
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In equations IV-22, Et-^, St^TD^, etc., refers to the sum of 
the time values, the sum of the products of time and time difference 
values, etc. , for the whole set of data available at the fixed point. 
Hie total number of sets of TD vs t data is n. 
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V. RESULTS AND OCNCLUSIQNS 



Part A of this chapter presents the results of calculations based on 
the theoretical principles previously discussed and data obtained from the 
experiment. Part B provides a summary of the results obtained, how these 
results differ from what would be expected based on the theory, possible 
reasons for any variations from expected outcome, and a general conclusion. 

A. Calculations 

1. Azimuths, Directions of Hyperbolas, and Directions of ID's 

Azimuths of the Loran-C stations from the survey points are one 
of the basic requirements of several calculations used in this 
experiment. Station positions of Range-7 and Luces Point (the primary 
stations used in this experiment) and other points were obtained from 
Horizontal Control Data published by the U.S. Dept, of Commerce, Coast 
and Geodetic Survey (see Appendix A) . Station positions of the Loran-C 
transmitters were obtained from Loran-C Data Sheets provided by the 
Defense Mapping Agency, Hydrographic/Topographic Center (see Appendix 
B) . For convenience in avoiding datum conversion, the North American 
Datum 1927 (NAD 27) was used because published positions based on this 
datum were available for all of the station positions used in the 
experiment. Station positions are summarized in Table V-l. 
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Azimuths ; Azimuths from Range-7 and Luces Point to Loran-C 
stations 9940-M, W, X, and Y were calculated using a Fortran 77 
subroutine provided by EMAOTC (Appendix C). (Operation of the, 
program was checked using kncwn data points and azimuths. Results of 
these computations are presented in tabular form in Table V-2. For 
convenience in visualizing the angular relation- ships, results are 
also presented in graphic form in Figures V-la and V-2 a. 

Directions of Hyperbolas : The directions of hyperbolas for 
9940-W, X, and Y were calculated using equations TV-4 and the azimuths 
obtained as described in the previous paragraph. A sample calculation 
follows: 



Hyperbola Direction for 9940-X at Range-7: 



HDr7 



(Azimuth of 9940-X) + (Azimuth of 9940-M) 



2 

343.071686 + 38.358551 
2 



190.715118° 

(from North) 



Decimal degrees are used in the calculations for convenience. 
Results of the calculations are again presented in Figures V-lb and 2b. 

Direction of Increase of TD : Directions of increasing TD are 
at right angles to the hyperbola directions and increase in the 
direction of hyperbolas which approach the master station. This 
involves either adding or subtracting 90" from the direction of the 
hyperbolas. For checking purposes, a Loran-C chart of the area was 
used. Results are again presented in Figures V-l and 2. 
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TABLE V-l. Station Positions (NAD 27) 



Station 


Latitude 


Longitude 


Range 7 


36° 39' 02.47787" N 


121° 49' 08.58202" W 


Luces Point 


36° 38' 10.524" N 


121° 55' 38.399" W 


9940-M 


39° 33' 07.046" N 


118° 49' 52.241" W 


994 0-W 


47° 03' 48.594" N 


119° 49' 52.241" W 


9940-X 


38° 46' 57.472" N 


122° 29' 40.050" W 


9940-Y 


35° 19' 18.342" N 


114° 48' 13.946" W 



TABLE V-2 . Azimuths Fran Range-7 and Luces Point. 



Range-7 to 


Degrees, min, sec. 


Decimal Degrees 


9940-M 


38° 21' 30.784" 


38.358551° 


994 0-W 


7° 32' 22.823" 


7.539673° 


9940-X 


343° 04' 18.070" 


343.071686° 


9940-Y 


101° 47' 36.229" 


101.793397° 


Luces Point to 






9940-M 


38° 51' 08.262" 


38.852295° 


994 0-W 


7° 45' 43.837" 


7.762177° 


9940-X 


344° 23' 49.981" 


344.397217° 


9940-Y 


101° 36' 16.319" 


101.604533° 
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9940-W 7.539673* 



9940-X 343.071686* 




9940-M 38.358551* 



9940-Y 101.793397' 



FIGURE V-la. Azimuths Fran Range-7 to 9940-44, W.X, and Y. 




FIGURE V-lb. Directions of Hyperbolas 



190.715118* 
at Range-7. 



9940-Y 




9940-W 112.949112* 

FIGURE V-lc. Directions of Maximum Increase (Gradient) of ID's at R-7. 
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9940-W 7.762177 



9940-X 344.397217* 




9940^4 38.852295* 



9940-Y 101.604533' 



FIGURE V-2a. Azimuths Fran luces Point to 9940-W,W.X,arc3 Y. 




FIGURE V-2b. Directions of Hyperbolas at luces Point. 



9940-Y 




9940-W 113.307236* 

FIGURE V-2c. Directions of Maximum Increase (Gradient) of ID's at IP. 
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2. Lanewidths 



Lanewidths were calculated from Equation IV-3a using the 
azimuths computed in Section V-A-l . Note that the units of lanewidth 
are m/iisec. With a standard lane of 1 nsec the width of a lane is the 
same number in meters. Velocity of transmission was calculated using c 
and the index of refraction specified in the table of Loran-C constants 

supplied by the Defense Mapping Agency (Ref. 4). Sample calculations of 
lanewidth follow: 



Velocity of Transmission = 



n 



2.99792453 x 10 8 m/sec 
1.000338 

m 

= 2.996911624 x 10 8 n^sec 

Lanewidth of 9940-Y at Range-7 for standard 1 n sec lane: 



l R7 “ 



(10 -6 sec) (2.996911624 X 10 8 m/s ec) 



101.793397* - 38.358551* 



2 sin 



[V-l] 



= 285.024 m 



Results of the calculation of lanewidths are given as the L 
values in Table V-3 . In the table, L^, Lq , and L 3 correspond to the 
lanewidths for 9940-W, X, and Y, respectively. 
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3. Average TD's, Standard Deviations, Variances 



Because of the large amount of data involved, it was necessary 
to compute statistical information using a computer program. Since 
data was in the form of three TD's (one for each master-secondary pair) 
every five seconds, and readings were taken for periods of 20 to 40 
minutes at each station, a Fortran 77 program with subroutines for 
statistical information was written (see Appendix C) . 

Station Range-7 was the "known" station of the differential 
system, at which more readings were taken than any other. For this 
reason, and because regression equations of TD vs Time were needed only 
at that station, statistical information was calculated only for 
station Range-7. 

Mean of the TD's was computed by the following (Ref. 18): 



ZTDp 

Mean TD - [V-2] 

n 



Variances of the TD's were calculated using the following (Ref. 18): 



ZTDj 2 - (STDp) 2 /N 

s 2 = [V-3] 

N - 1 



Standard deviations were calculated from variances by (Ref. 18): 



formula (Ref. 18): 

s = (s 2 )^ [v-4] 
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Equations V-2 , V-3, and V-4 were incorporated into a Fortran 77 
computer program which was used to calculate statistics of TD's for 
9940-W, X, and Y. 

It should be noted that these statistics were in terms of time 
differences, while in same cases it was necessary to knew standard 
deviations and variances in terms of distances. Equations IV- 5 were 
used for this conversion where necessary. 

Corrected vs Uncorrected Data : One use of the statistics of 
the data was to compare the standard deviation of raw data to the 
standard deviation of data which was corrected for three secondary 
stations. It was considered unlikely that the three station correction 
could remove bias from the TD readings (error in terms of absolute 
position) because most of this is due to systematic error rather than 
random variations. It was believed that the three station correction 
would be most effective in reducing the effect of random variations in 
time differences. The statistics on uncorrected and corrected data show 
results consistent with this expectation (Table V- 7). 



4. Error Ellipse and Circular Error at Range-7 

Range-7 was used as the known point of the differential system, 
and a larger amount of data was collected at that point. Using the 
crossing angles, standard deviations and variances of 9940-X and Y (the 
best pair at that point) , dimensions and orientation of the standard 
error ellipse and CEP 50% and 90% circle radii were calculated for the 
point using the information outlined in section III-B. 
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The standard error ellipse and error circles for 50% and 90% 
were plotted over a scatter plot of data from Range-7. The result is 
shewn in Figure V-3. 

Calculations of the standard error ellipse dimensions were as 

follows: 



Using 9940-X,Y: 

°1 = a 9940-Y = (0.072546) (285.024) = 20.677m 

O 2 = a 9940-X = (0.068638) (322.024) = 22.103 m 

0 = 250.075974° - 190.715118° = 59.360856° 



Calculating the semi-major and semi-minor axes of the ellipse: 



2 



20 . 677 2 +22 . 103 2 + ( (20 . 677 2 +22 . 103 2 ) 2 -4 (20 . 677 22 2 . lO3 2 sin 2 0) ) ^ 



■*x 



■'x 



2 



2sin 2 0 

= 30.595 m 

20 . 677 2 +22 . 103 2 - ( (20. 677 2 +22 . 103 2 ) 2 -4 (20. 677 22 2 . lO3 2 sin 2 0 ) ) % 



2sin 2 0 



= 17.362 m 



Orientation of the major axis of the ellipse: 



tan 2© = 



20.677 2 sin 2(59.360856°) 



20.677 2 cos 2(59.360856°) + 22.103 2 
© = 26.474268° 



Azimuth of major axis = 250.075974° - 26.474268 = 223.601716° 

= 223.601716° - 180° = 43.601716° (from North) 
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1 100 
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i 



100 ^ 



scale in meters 



i 



Figure V-3. Scatter Diagram, Standard Error Ellipse, and CEP^q and 

CEPg 0 Error Circles, for Range-7 Data (see text for details) . 
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Using tables in Reference 6 , pp. 18-20, radii of the 50% and 
90% error circles were calculated as: 

CEP 50 = 27.936 m CEPgo = 63.226 m 

Because the data included in the plot of Figure V-3 includes 
two intervals of data collection, the points are not as evenly 
distributed as they would otherwise be. For this reason there are 
fewer points near the center of the plot than would normally be the 
case. The error circles and ellipse seem to be reasonably correct for 
the given data. 

5. Correction Factor for Three Secondary Stations 

The correction for three secondary stations was designed as a 
technique to reduce the random error associated with time difference 
data. It was expected to increase the precision of the position rather 
than the absolute position. 

For a variety of reasons there will be some covariance between 
concurrent time difference readings from different master-secondary 
pairs. The correction for three secondary stations is not expected to 
reduce the error associated with this type of randan error. The non- 
correlated randan error in TD readings should be reduced by this 
correction to sane extent. 

In order to reduce the number of computations and the time 
required to process data, most constant factors associated with the 
correction factor can be computed in advance and used as constants in 
the computer program. This reduces the correction to only a few lines 
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of relatively simple computation. Because of the large amount of data, 
this reduces significantly the time required to compute corrections. 

In the calculations of adjustments to positions of Range-7 and 
Duces Point, the only factors which differ significantly are the long 
term time differences at the two points. In this case the same A_^ 
factors can be used for correcting both positions for three seco ndar y 
stations. In cases where seme of the parameters may differ 
significantly, it would be necessary to calculate different A^ factors 
for each position. Table V-3 gives parameters for input to the Fortran 
77 correction program for calculation of the A-[ parameters and 
calculation of the corrections to data. 

TABLE V-3. Parameters used in calculation of correction 

factor k . 



Parameter 


Value 


*1 


112.949112° 


0 2 


100.715118° 


0 3 


340.075974° 


CT 1 


0.150606 0 sec 


a 2 


0.068638 /xsec 


a 3 


0.072546 /xsec 


Ll 


563.934 m//isec 


^2 


322.966 it v//isec 


*3 


285.024 m//isec 
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6. Regression Equations for Time Difference vs Time 



Since only one Loran-C monitor receiver was available, it was 
necessary to develop a method of estimating TD's at the known point 
during times when differential corrections were required at the unknown 
points. Linear regression equations of TD vs time were used for this 
purpose. 

Using equations IV-21 and IV-22, the data from Range-7 was 
processed to obtain the required summations of time and time difference 
necessary to calculate the b 0 and b^ coefficients of the linear 
regression equation. Only Range-7 data from the nearest time periods 
before and after the differential correction period for Luces Point 
data were used. 

The procedure for processing the data was first to apply the 
correction for three secondary stations to the Range-7 and Luces Point 
data during the time periods of interest. Range-7 data was then 

processed to obtain the necessary summations for the calculation of the 
b parameters of the regression equations. Using the regression 
equations for Range-7 data and long term TD averages at the two points, 
differential corrections for the Luces Point data were calculated as 
described (Section V- 7). Fortran subroutines (Appendix C) were added 
added to the statistics program to obtain the summations and a separate 
program was used for the differential correction of Luces Point data. 
Results of the corrections were compared to calculated values of TD's 
provided by the Defense Mapping Agency, Hydrographic/Topographic 
Center. TD vs time is plotted in Figures V-4 along with the regression 
line for the data. Regression equations are in Table V-4. 
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TIME vs TD 

9940-W 9/28/07 



16316.1 .. 




16316.1 



+ 



9.Q 



__l 1 1 1 h 

10. a 11.0 13.0 13.0 14.0 

TIME (PST) 



Figure V-4a. TD vs time and the regression line for 9940-W 
for Range-7 data corrected for three secondary stations. 
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Figure V-4 b. TD vs time and the regression line for 9940-X 
for Range-7 data corrected for three secondary stations. 
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Figure V 4c. TD vs time and the regression line for 9940-Y 
for Range-7 data corrected for three secondary stations. 



63 



TABLE V-4 . Regression Equations of TD vs Time at Range-7. 



994 0-W: 


TDw = 


16316.51564 


- 0. 00807243t 


9940-X: 


TDx = 


27522.33725 


- 0. 01385199t 


994 0-Y: 


TDy = 


42749.29371 


+ 0. 010567 07t 



time units: hours with minutes and seconds in decimal hours. 

7 . Differential Corrections 

Eguation IV-20 describes the differential correction which is 
applied to TD readings at an unknown point based on TD's at a known 
point. In this experiment the positions of all points are known, based 
on their published positions from survey data (see appendices) . One 
point, Range-7, was chosen as the known point for purposes of 
differential corrections. Data from one other point was chosen as the 
unknown point for purposes of comparison to its known position. 

Three things must be known to calculate a differential 
correction based on eguation IV-20. Rather than use long term averages 
in this eguation, the expected TD's at both points were computed by a 
standard Loran-C program. Current TD's at Range-7 were calculated from 
the regression equations described in the last section. The ASF 
corrections in equation IV-20 were calculated for each point by the 
Defense Mapping Agency using a computer program which is used in the 
preparation of ASF correction tables, but which can be used for 
determining the correction at individual points. The difference in ASF 
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corrections is used in equation IV-20. Table V-5 gives expected TD 
values (used as long term averages) at the two points of interest: 

TABLE V-5. Calculated TD's used as long term averages (in nsec) . 





Range-7 


Luces Point 


9940— W 


16314.68 


16300.09 


994 0-X 


27523.13 


27496.79 


994 0-Y 


42749.49 


42755.87 



ASF corrections were approximately the same at both points due 
to the nearness of the points. Differences in ASF corrections used in 
equation IV-20 were very small: 0.005, 0.002, and 0.003 nsec for 9940- 
W, X, and Y, respectively. Corrections of this magnitude are probably 
less than the accuracy of the data used to calculate them, but at 
longer ranges these may be large enough to be significant. 

The final equation for the differential corrections based on 
TD's at Range-7 are in Table V-6. Note that t is Pacific Standard Time 
in hours, with minutes and seconds expressed as a decimal. 

TABLE V-6. Differential Corrections (C^ in jxsec, t in hours PST) 



< -dw - 


-1.83064 + 


0. 00807243t 


<^x “ 


0.79575 + 


0.01385199t 


^dy ~ 


0.19829 - 


0. 01056707t 
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8. Accuracy of Position Fix 



Figure V-5a shews a plot of uncorrected Range-7 data based on 
expected TD's calculated frran the known geographic position. The grid 
is in meters, oriented north-south and east-west, and with the center 
at the published geographic position of Range-7. When the Correction 
For Three Secondary Stations is applied and the average of the Range-7 
data is placed at the origin, the result is as shewn in Figure V-5b. 

Note that the plot is new centered on the origin due to the 
fact that all absolute error (bias) has been removed from the data. 
The three station correction has also decreased the standard deviation 
and variance of the Range-7 data as shewn in Table V-7. 

Figure V-6a shows Duces Point data before application of any 
corrections. These points are plotted based on the expected TD's at 
the knewn coordinates of the point. The center of the uncorrected 
Duces Point data is at x = 375.820 m, y = 85.501 m, and at a range of 
385.423 m and at an azimuth of 77.183° (from North) from the origin. 
Radii of the CEPgQ and CEPgQ error circles are 14.9 m and 28.4 m, 
respectively . 

Figure V-6b shows Luces Point data after the application of 
both three station and differential corrections. Again, the origin is 
at the knewn position of the point. For visual comparison purposes, 
the same scale is used for both Figures V-6a and b. For the corrected 
data, the center of the data is at x = 47.213 m and y = -7.015 m, with 
a range of 47.731 m and at an azimuth of 98.451° (from North) from the 
origin. CEP 50 and CEPgg radii were 12.6 m and 26.0 m, respectively. 



66 




Figure V-5a. Scatter Diagram For Range-7 Uncorrected Data. Scale is 
in meters oriented N-S and E-W, centered on the time 
differences calculated from the published coordinates 
of the point. 
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CORRECTED 9940-X, Y 
SCATTER DIAGRAM 
STATION RANGE-7 
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Figure V-5b. Scatter Diagram for Range-7 Data Corrected for Three 
Secondary Stations and Differential. Scale is in 
meters oriented N-S and E-W. Center of data is 
corrected to correspond with calculated time 
differences for the published coordinates of the point. 
See text for details. 
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Figure V-6a. Scatter Diagram For Luces Point Uncorrected Data. Scale 
is meters oriented N-S and E-W, centered on the time 
differences calculated from the published coordinates 
of the point. 
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Figure V-6b. Scatter Diagram for Luces Point Data With Differential 
Corrections Based on tenge-7 Data and With Correction 
For Three Secondary Stations. Origin of plot is at the 
time differences calculated from the published 
coordinates of the position. 
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TABLE V— 7. Standard Deviation and Variance of Range-7 Data With and 
Without the Correction For Three Secondary Stations. 



Pair 


Without Correction 


With Correction 


o 


0.068638 


0.051161 


9940-X 






a 2 


0.004711 


0.002617 


a 


0.072546 


0.052081 


9940-Y 






a 2 


0.005263 


0.002712 



B. Conclusions 

Based on the results given in part A, the absolute accuracy of 
the Duces Point Data was improved from an average error of 385.423 
meters to an average error of 47.731 meters. Practically all of this 
improvement was a result of application of the differential correction. 

Precision of the Duces Point data was improved from 14.9 m and 
28.4 m for the 50% and 90% CEP radii, uncorrected, to 12.6 m and 26.0 m 
for the corresponding radii, corrected. Most of this improvement is a 
result of the correction for three secondary stations. 

These figures are based on the processing of 486 sets of TD 
data from Range-7, the "known" point of the differential system, and 
252 sets of TD data at Duces Point. Each set of data included time and 
three TD's, one for each master-secondary pair (9940-W, X, and Y) . 

Under ideal conditions a better absolute correction would be 
expected for these points. From an examination of the characteristics 
of the data, the weakest point is probably the lack of a continuous set 
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of TD readings at Range-7. Error in the estimate of TD's based on the 
linear regression lines may be as much, or more than, the 48 -meter 
average error in the corrected positions calculated for Luces Point. 

Other sources of error may include the fact that seme of the 
Loran-C signals traversed seme seawater path near Luces Point, with 
unknown results on TD readings at that point. This type of error is 
known to result under similar conditions, but no means of estimating 
its magnitude is available at this time. 

In summary, the absolute positions calculated from differential 
corrections of Loran-C data seem to be a significant improvement over 
the positions calculated from raw data. While the correction for three 
secondary stations does not add much to the estimate of positions in 
this experiment, the improvement of 50% CEP radii by about 15% does 
seem significant enough to warrant further consideration. 
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APPENDIX C-l 



C THIS PRO&RAH READS OATA FROH A FILE OF TIHE, A, 6, C 

C CONTINUES TO THE END OF THE FILE, CALCULATING THE HEAN, 

c variance, and standard oeviation of the a, b, AND C COLUflNS, 

C AND PRINTS OUT THE RESULTS INTO A FILE tdout.f 

double precision t i ae , t i ae 1 , t i ae2, t i ae 3 , e , b , c , na , na , nb , nc , 

4sua, suae, suab, suae. su#e2,sueb2,suec2.suael.$uabl,suacl,al,bl, cl 
open (unit«l, file * ’out9.f’, for* * ’foraetted', 

4access=’ sequential’) 

C open (unit*2,file * 'outfc.f', fora * ’foraetted 1 , 

C ♦access* ’ sequent i a 1 ’ ) 

rewind 1 

C reed! 1 , ■ ) f rst In 

C re a d ( 1 , * ) secln 

SO write (*.4) t i ae 1 , t i ae 2, 1 1 ae 3 , a 1 , bl , c 1 , na , na , nb , nc 
tiae » tiael * tiae2/60. ♦ tiae3/360Q. 

C write (2,8) t i ae , a 1 , b 1 , c 1 

tuft * tua 4 1. 
sura c suae 4 a 
sural a sural 4 el 

$ure2 * sure2 ♦ a*e 

suab * SUBb 4 b 
surbl a surtl 4 bl 

sueb2 a suab2 ♦ b*b 

sure a sure ♦ c 
suae 1 a suae 1 4 Cl 

SUt.c2 a SU*C2 4 c*c 

rtad( 1 ,2, errxio. iostet*ioerr) tiael,tire2,ti»e3,el,a.bl,b, 

♦ c 1 , c 

if I loe-r) 10.50,10 
10 sur = suit - 1. 
wr i t e ( * . 6 ) ioerr 
call s 2 ( s u a a 2 . sure, sur .sural) 
call s?(surb2, suab, sur, surbl) 
call s2 ( sumc2 . suae , sur , suae 1 ) 

2 for rat It 1 , d 2 . 0 , 1 4 , O*.0,t7,d2.C,tl£,d£. 3 , t 2 1 , d 5 . 3 , 1 2 7 , 
4d8.3,t30.d5.3,t36.d£-3.t39,d5.3) 

4 foraat(f3.0,2*,f3.0,2»,f3.0,2i,fl0.3,2»,fl0.3,2*,fl0.3,2*, 

4iA, 1 4 , iA, 14) 
t for«Bt ( ’ ioerr a * , i 3 ) 

8 for»at(fl2.£.2»,fl0.3,2i,fl0.3,2i,fl0.3) 

9 for»at(‘deta error’) 

40 stop 

end 

C 

subroutine s 2 ( i 2 , i , d , s 1 ) 

double precision t2, i ,d, aean, varian, sd, si 
a e an a sl/d 

varian a ( i 2 - i * i / d ) / { d - 1 ) 
sd * dsqrt(verian) 
write(*,12) ae an , v ar i a n , sd 
C wr i t e ( 2 , 1 2 ) ae an , v a r i an , sd 

return 

12 f or »a t ( ’ aean « * , f 1 2 . 3 , 2 r , ’ ver i a nee a ’ , f l 2 . 8 , 2 * , ’ sd » ' .fl2.fi) 

end 
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10 

20 

50 

AO 

50 

60 

70 

80 

90 

95 

100 

105 

no 

120 

130 

1 A 0 
150 
160 
170 
175 
180 
185 
190 
200 
210 
215 
220 
22 1 

2 3C 
240 

'250 
255 
260 
2t 5 
270 
280 
290 



re • 

re* --- PROGR Art FOP 6 AUS SIAM El 1 rt IN AT 1 OK U1TH PIVOTAL CONDENSATION 

rei FOP SOLUTION OF SYSTErtS OF LINEAR EQUATIONS 

res U1TH UP TO 10 UNKNOWNS 

re* 

re* 

open 'i*. II, ’ytest.f’ : di* 
input II, rows ! , cols! 



cols! 



input the atstrii: 



r ■ rows! 
re* 

re* 

re* 

for j ■ 1 to rows ! 
for k « 1 to cols! 
input 1 1 , A » ( j , k ) 
if k <> cols! then 160 
A I ( j , k ) - >1 1 All j.k) 
ne 1 1 k 
ne 1 1 j 
re* 

re* ----- 
re* 

gosub 500 
gosub 890 
gosub 650 
re* 

if h l g h r o w 
re nr 

gosub 700 : 
g o s u t 7 1 0 : 
gosub 820 : 
r e xi. 

rowl = rowl 
ret 

if rowl (> rows then 200 : 
gosub 760 : re* -- 

gosub 890 : re* -- 



AI(10, 1 lhdi* B • ( 1 1 ) : d i ■ ni(n):di* ii(ll) 
» 1 : col 1 « 1 : highrow « l 



begin *ein program 

re* 

re* 

re* 

rowl then 2A0 

r e i. 

r e a 

r e » 



check if uniaue solution 
prints out array 

checks for row with highest first non-zero entr 



eichenge rows if necessary 

sets first non zero entry in rowl to 1.0 

sets first entry to zero in following rows 



coll * coll 



highrow s rowl 



r eir 



-- increment count f 



rert --- checks if lest row has beer, reached 
sets first non zero entry in final row to 10 
final printout of array 



295 re* 

300 ren calculates the unknown x velues: 

305 re* 

310 i • ( c o 1 s ! ) * 1 : a * rows! 

320 for i ■ a to 1 step -1 
330 for j s i to cols! - 1 
3 4 0 I I f i ) r i | ( j ) - I I ( j ♦ 1) * A I ( i , j ♦ 1 ) 

350 next j 
360 ne i t i 
365 re* 

370 re* prints out the i values: 



375 ret 

380 for k * 1 to rows! 

390 print * B ' k-1 * * * ilfk) 

A 0 0 n e 1 1 k 
A10 end 

A20 re* end of *ain program. 

A 2 5 re* 

A30 re* begin subroutines - 

A 3 5 re* 

500 for i * 1 to rows! 

510 for j * 1 to rows! 

520 if i * j then 610 
530 for k » 1 to cols! 

5 A 0 if A • ( i , k ) « 0 t hen 560 
550 HI Ik ) * A»(J,k)/Al(l,k) 

560 ne i t k 
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570 for i > 1 to col*! - 1 

580 if ni(a) <) M(»*l) then 610 

590 neit • 

600 print ‘E0UAT10N5 NOT IN0EPEN0ENT ... NO UNiOUE SOLUTION ... ONE SOLUTION F0LL0U5* 

610 ne i t J 

620 neit i 

630 return 

635 re* 

640 rei ----- subroutine to check for roy with highest first non-tero entry: 

650 for i * rowl to rows! - i 

660 if ibsUKi.coll )) ) ebs ( Al ( i ♦ 1 , col 1 ) ) then 680 
670 highrow ■ i ♦ 1 
680 neit i 
690 return 

695 rei »unroutine to eichonge rows end put highest first entry first: 

700 f or i * col 1 to cols ! 

710 BK i) * AKrowi , i ) 

720 AKrowi, i) * Al ( hi ghrow , i ) 

730 A I ( h i ghrow , i) * Bt ( i ) 

74 0 ne 1 1 1 
750 return 

755 rea subroutine to set first entry in first row to 1.0; 

760 FI * AKrowi , col i ) 

770 for k * coll to cols! 

780 if AKrowi, coli) * 0 then 310 

790 AKrowi. k) * Alfrowl, k)/RI 

791 print k 
800 neit k 
810 return 

815 rei subroutine to set other first entries to zero: 

82C for j * rowl ♦ 1 to rows! 

830 PI * At f j, coll )/Al(rowl, coll) 

84 C for k s coll to cols! 

850 A«(J.k) * Al(j.k) - Alfrowl.k )*PI 
860 neit k 
870 neit j 
880 return 

685 res subroutine to print out the aatrii A; 

890 for J « 1 to rows! 

900 for k * 1 to cols! 

910 print A t ( J , k ) , 

920 neit k 
930 print 
940 neit j 

950 print rowl; * - - - 

960 return 

970 print ‘NOT ENOUGH E0UAT10NS ... HANY SOLUTIONS ARE POSSIBLE* 

980 return 
990 end 
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PROGRAM A 'n 
A * t 372206 4 

B r 6256563.6 

OPENfUNIT*] . riLE='o:in. f ' ,FORn=’FORnAT7EO* , 

*ACCESS=’SE0UENT1AL' ) 

REWIND 1 

C erimutN is fror point 1 to point 2 

readi i.2) alat:,alon: 

READfl , 2 ) Al ATI , A l ON 1 
FLAG * 0.0 

CALL lNVlA.B.ALAn.AL0Ni,ALAT2,AL0N2.FLAT,DlST.AZ) 

A? s AZ ♦ 3. 141S<>?654 
AZ = AZ/3. 141592654*160. 0 
IF I A ? - 360 . 0 ) 10,10.20 
20 A 2 * A Z - 36 C . 0 
10 WRITE!* ,3) DIST.AZ 

X = ALAT 2/3. 1 4 1 592t 5<. *16 0.0 

Y * ALONZ/3 . 14 1 5*2t$4 1 180.0 

V = ALAT l /3 141592*54*160 0 
U = ALOM/3. 141592654*160.0 
WRITE » * . 4 ) * . Y 

WRITE f * . 5) V. L 

2 FOFHAT f F<5. E . T 1 2. F9. £ I 

3 FORP.aT f ’DISTANCE = * . F 1 2 . 3, 1 * . ’T1ETER5* . 3> . ’ AZinL'TH = * . F 1 2 E . 2 X . 

4* | NORTH 1 ’ ) 

4 FORr.ATl’FCJNT;: ',2 x,‘LAT:*, 2X.FE.5, ‘N*,3x.’LON.’,2>.F9.5. .’W*1 

5 FORMAT {’ F01 NT 1 : ’ . 2 > , * L a t . * , 2 x . F l 5 . * N * . 3 > . ’ L OK : ’ . 2 » . F c 5.*W) 

E NO 

C 

suer Du tine IN.MEORAI. fofad.rlti . rlm .plt:.rln:.cflc*.di 5T. A21 
I c I C ' L G 1 101. 1C 1.102 
101 FLA’ = l O-FORAC 'El RAD 
Fa. at: - F L A T ’ r L A 7 

FJ = Falat:* 1 25 

f : = f a L a • : • c . 5 
F 3 = FAL AT2*C. :* 

F 4 S F A L A T 2 * 0 . 1 : 5 
F 5 - FALAT 2*0. 0 1 : r 
Ft s FALAT2 + F1.A 1 
F - = F t * 1 . 0 
F6 = Ft *0. 5 
F 1 = 3.141 592654 

TWPFI = 6. 263165 3 CT 
C F L G * 1.0 

1 0 : EEIA1 = ATAn((1.0-FLAT)’SIN(RlT1)/C05IRLT1)) 

SEE T a i = SIMEETAl) 

CEETAl = COS f E E " a l ) 

BETA! = ATAf.f (l.0-FLAT)*SlN(RLT2)/C0SrRLTD) 

SEETA2 = S 1 N ( B E T A : ) 
ceeta: s c o s i e e t a : } 

DELL s FLN1-RLN2 

AOEIL = AESICELLl 

IF 1 ADELl-FI ) lOt.lC5.lG' 

105 ACELL = TUOFI-ACELL 
1C* S1DEL = SIN.ACElLI 
CODEL = COS! ADEL L) 
a t SEE T a i • SEE T A2 
e = cbeta: *ceeta: 

COFHI = A 4 6 * C ODE L 
SIFHI = 

8S0FT [ f SIDEt *CEE TAT 1 * ' 24 [ SBETA2*CEETA1 -SEETA1 'CEETA: * CODEL I * ’ 2 i 
C = B*S1DEL/SIFHI 
Er = 1 . C-C * C 
Fhl = A S I N I S 1 P M I 1 
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IF (COPMI) 207,206.20ft 

207 PHI * PI-PHI 

208 PHI 50 * PHI * PH 1 
CSPHI « 1.0/SIPH1 

C T PHI * C0PHI/5IPHI 
PSYCO * SI PHI * COPHI 
TERH1 * F 7 * P H I 

TERH2 * A*(Ffc*SIPHl-F2*PHlSQ*C$PHI ) 

TERH3 * En* (F2*PHlSQ*CTPHl-FB* (PHI *PS YCO ) ) 

TERH4 » A*A*F2*PSYC0 

TERMS * En*EH*(FS* (PHl4PSYC0)-F2*PHISQ*CTPHl-F4*PSYC0*C0PHl*C0PHl ) 
TERH6 - A*fH*F2*(PHI$Q*CSPHl4PSYCO , COPHl) 

0 1 ST * P0RADMTERni4TERn2«TERH3-TERH4*TERnS«TERH6) 

T E RH 7 * F 6 * PH I 

T E RFlfi * A* (F2*$IPHl4FALAT2*PHlSQ*CSPHl ) 

TERH9 « EH* (F3 , F3*PSYC0*FALAT2*PH1SQ , CTPHI-F1*PHI ) 

ZLAH - C* (TERH7-TERH8*TERH9)*A0ELL 

CTAZs(5BETA2*CBETA1-C0S(ZLAHZ) *5BETAI *CBETA2)/(S1N(ZLAH) ‘CBETA2) 

IF (CTAZ) 210.209,210 

209 CTAZ * 0.00000005 

210 AZ * A T AN ( 1.0/CTAZ) 

IF 1 DELL ) 215,214,214 

214 IF (DELL-PI) 211.212,212 

211 IF (CTAZ) 220,221,221 

220 AZ * A Z ♦ P I 
GOTO 221 

212 IF (CTAZ) 217.2ie.21C 

2 1 5 IF ( DE L L ♦ P I ) 211,211.216 

216 IF (CTAZ) 217,218.218 

217 AZ * PI-AZ 
GOTO 221 

218 AZ = TUOPI-AZ 

221 AZ * AZ 4 F I 

AZ * AZ-TUOPI 
IF (AZ) 213,219, 219 

213 AZ = AZ+TU0P1 

219 RETURN 
END 
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